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The effects of perturbative Lorentz and CPT violation on neutrino oscillations are studied. Fea- 
tures include neutrino-antineutrino oscillations, direction dependence, and unconventional energy 
behavior. Leading-order corrections arising from renormalizable operators are derived in the general 
three-flavor effective field theory. The results are applied to neutrino-beam experiments with long 
baselines, which offer excellent sensitivity to the accompanying effects. Key signatures of Lorentz 
and CPT violation using neutrino beams include sidereal variations in the oscillation probabilities 
arising from the breakdown of rotational symmetry, and CPT asymmetries comparing neutrino and 
antineutrino modes. Attainable sensitivities to coefficients for Lorentz violation are estimated for 
several existing and future experiments. 



I. INTRODUCTION 

Experimental investigations of neutrino properties 
have provided crucial insights into particle physics since 
the existence of neutrinos was first proposed in 1930 by 
Pauli [1] to explain the spectrum of beta decay. In re- 
cent years, the confirmed observation of neutrino oscilla- 
tions has established the existence of physics beyond the 
minimal Standard Model (SM) [2]. The interferometric 
nature of the oscillations makes them highly sensitive to 
new physics, including potential low-energy signals that 
may originate in a fundamental theory unifying quantum 
physics and gravity at the Planck scale mp ~ 10^^ GeV. 

In this work, we investigate the experimental implica- 
tions for neutrino oscillations of Lorentz and CPT viola- 
tion, which is a promising category of Planck-scale signals 
[3] . The SM is known to provide a successful description 
of observed phenomena at energies well below mp. As 
a consequence, the manifestation of Planck-scale effects 
involving Lorentz and CPT violation is expected to be 
well described at accessible energy scales by an effective 
field theory containing the SM [4, 5]. 

The comprehensive effective field theory describing 
general Lorentz violation at attainable energies is the 
Standard-Model Extension (SME) [6, 7]. It incorporates 
both the SM and General Relativity, serving as a real- 
istic theory for analyses of experimental data. In the 
SME Lagrange density, each Lorentz- violating term is an 
observer scalar density constructed as the product of a 
Lorentz- violating operator with a controlling coefficient. 
Under mild assumptions, CPT violation in effective field 
theory is accompanied by Lorentz violation [8], so the 
SME also describes general breaking of CPT symmetry. 
These ideas have triggered a wide variety of tests over 
the past decade [9]. Several experimental searches have 
been performed using neutrino oscillations, yielding high 
sensitivities to SME coefficients for Lorentz and CPT vi- 
olation [10-12]. 

Since both Lorentz- violating operators and mass terms 
can induce neutrino mixing, one way to classify neu- 
trino models with Lorentz and CPT violation is in terms 
of their neutrino-mass content. Three categories exist: 



massless Lorentz- violating models, in which no neutrinos 
have mass; hybrid Lorentz-violating models, with mass 
terms for a subset of neutrinos; and massive Lorentz- 
violating models, where all neutrinos have conventional 
masses. 

In massless Lorentz-violating models, all observed neu- 
trino oscillations are attributed to nonzero coefficients for 
Lorentz violation rather than to masses. Certain coeffi- 
cients can combine via a Lorcntz-seesaw mechanism to 
produce pseudomasses that mimic the behavior of mass 
terms for a range of neutrino energies [13]. The proto- 
typical example is the bicycle model [14] which uses two 
nonzero coefficients for Lorentz violation to reproduce 
the expected behavior of atmospheric neutrinos. This 
model agrees well with atmospheric data from the Super- 
Kamiokande experiment [11]. However, a combined anal- 
ysis of neutrino data excludes both the bicycle model in 
its simplest form and a five-parameter generalization [15]. 
Massless models may also predict sidereal signals arising 
from the violation of rotation invariance [16]. The corre- 
sponding coefficients for Lorentz violation have been con- 
strained in experimental analyses by the Liquid Scintilla- 
tor Neutrino Detector (LSND) [10] and the Main Injector 
Neutrino Oscillation Search (MINOS) [12]. At present, 
it is an interesting open challenge to construct a massless 
Lorentz-violating model that is globally compatible with 
existing neutrino data. 

For hybrid Lorentz-violating models, the tandem 
model [17] is the sole exemplar. It is a three-parameter 
model containing one neutrino mass, one coefficient for 
CPT-even Lorentz violation, and one coefficient for CPT- 
odd Lorentz violation. The model appears globally com- 
patible with existing experimental data, including the 
LSND anomaly [18]. The tandem model predicted a low- 
energy excess in the Mini Booster Neutrino Experiment 
(MiniBooNE) prior to its discovery, although the ob- 
served excess is quantitatively greater [19]. This success 
suggests that further theoretical investigations of hybrid 
Lorentz-violating models would be of definite interest. 

Massive Lorentz-violating models are the primary fo- 
cus of the present work. Most existing data from neu- 
trino oscillations are consistent with oscillation phases 
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proportional to the baseline L and inversely proportional 
to the energy E. This is conventionally interpreted as a 
consequence of mixing induced by a nondegencrate mass 
matrix. In massive Lorentz-violating models, the mix- 
ing due to mass is assumed to dominate over that due 
to Lorentz violation. Our goal here is to present a gen- 
eral study of perturbative Lorentz and CPT violation on 
mass-induced mixing, valid over a wide range of L and 
E. 

The analysis presented here incorporates all coeffi- 
cients for Lorentz violation associated with quadratic op- 
erators of rcnormalizable dimension in the neutrino sec- 
tor [13]. Using notation reviewed in Sec. II, these coef- 
ficients are {aL)tb^ {cl)^^ , S^f, and H^j^. Both {aL)'^^, 
and q"^ also control CPT violation. Taken alone, the 
coefficients (ar)";, and il"^ generate oscillation phases 
proportional to L but independent of E, while (cl)'^^ 
and g^^ produce phases proportional to the product LE. 
This indicates that experiments with long baselines or 
high energies are of special interest for studies of mas- 
sive Lorentz-violating models. However, the techniques 
outlined in this work apply for any baseline for which 
the perturbative approximation is valid, including ones 
where oscillations due to mass are negligible. Indeed, 
the expressions for oscillation probabilities presented here 
reduce to those obtained for massless Lorentz-violating 
models [16] in the limit of vanishing mass mixing. 

Massive Lorentz-violating models can exhibit effects 
lying in any of the six classes of physical effects due to 
Lorentz and CPT violation [13]. All coefficients affect the 
spectral dependence in at least some part of the energy 
range. Many of the associated operators violate rota- 
tion invariance, which can produce direction-dependent 
oscillations. For some experiments, including ones with 
neutrino beams, the daily rotation of the Earth induces 
variations in time of the probabilities at multiples of the 
sidereal frequency. Both CPT violation and neutrino- 
antineutrino mixings can occur. 

In this work, we show that for massive Lorentz- 
violating models the coefficients {aL)ab and (cl)"!! 
primarily affect neutrino-neutrino and antineutrino- 
antineutrino mixings, with (ai,)"^ controlling first-order 
differences between the two mixings due to perturba- 
tive CPT violation. Since the original introduction of 
these SME coefficients [6] , a substantial theoretical liter- 
ature has developed concerning their implication for neu- 
trino behavior in the context of massive Lorentz-violating 
models. Many works restrict attention to the special 
isotropic limits with only (01,)^^ i'^i^)ab nonzero and 
real [20-33], and in some cases also to two flavors. A few 
consider also anisotropic effects [13, 34-38]. Here, we 
treat the general case, allowing all components of (aL)^^ 
and {chYah *° ^® nonzero. 

In contrast, the dominant effects from g"^ and H% 

. . . ab ab 

in massive Lorentz-violating models arise only at second 
order. They involve neutrino-antineutrino mixing and 



also nonconscrvation of lepton number. A single flavor 
can therefore suffice to produce oscillations. Indeed, a 
simple analytical form is known for the mixing proba- 
bility for the general one-flavor case including mass [13]. 
A few theoretical works have considered special massive 
Lorentz-violating models of this type [39, 40]. At present, 
there are no published experimental constraints on any 
of the coefficients q"^ and H%. In this work, we investi- 

•^ab ab ' 

gate the general case and identify some potential signals 
for experimental searches. 

Overall, most massive Lorentz-violating models remain 
viable. Only a few percent of the available coefficient 
space has been explored experimentally [10-12]. The 
methods described in the present work demonstrate that 
access to essentially the whole coefficient space is avail- 
able via a combination of existing and future experi- 
ments. 

This paper is organized as follows. The basic theory 
and notation is presented in Sec. II. Section II A reviews 
the properties of the hamiltonian governing Lorentz and 
CPT violation in neutrino oscillations. The perturba- 
tion series for the transition amplitude is derived in Sec. 
II B, while the resulting oscillation probabilities arc ob- 
tained in Sec. II C. Section III considers ffist-order effects 
involving the coefficients (ai,)" and (cl)"'^- The direc- 
tional and sidereal dependences of the probabilities are 
discussed in Sec. HI A. Examples are provided for the 
case of three generations and its two-generation limit in 
Sec. IIIB. Asymmetries characterizing violations of the 
discrete symmetries CP and CPT are discussed in Sec. 
Ill C. Section IV investigates the second-order effects in- 
volving the coefficients g""^ and H". Oscillations with 
lepton-number violation are studied in Sec. IVA, while 
others arc considered in Sec. IV B. Section V concludes 
with a summary. 

II. BASIC THEORY 

This section begins with a brief review of the descrip- 
tion of Lorentz and CPT violation in neutrino oscilla- 
tions, assuming three generations of left-handed neutri- 
nos and their antineutrinos. We then use time-dependent 
perturbation theory to derive expressions for the tran- 
sition amplitudes and oscillation probabilities valid for 
small Lorentz and CPT violation. 



A. Hamiltonian 

Violations of Lorentz and CPT invariance in oscilla- 
tions of left-handed neutrinos and their antineutrinos can 
be characterized by a 6x6 effective hamiltonian {heff)AB 
taking the form [13] 

{heff)AB = {ho)AB + SHab- (1) 

Here, describes conventional Lorentz-invariant neu- 
trino oscillations, while 6h includes the Lorentz-violating 
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contributions. The uppercase indices take six values, 
A, B, . . . = e, 11, T, e, p., f, spanning the three flavors of 
neutrinos and antineutrinos. 

Under typical assumptions, the conventional term ho 
induces no mixing between neutrinos and antineutrinos. 
It is therefore block diagonal, and we write it as 



{h(-))ab 







1 / AmJ, 







(2) 



where E is the neutrino energy, lowercase indices 
a,b, . . . — e, indicate neutrinos, and lowercase barred 
indices a, &, . . . = e,p,f indicate antineutrinos. The two 
3x3 mass matrices are related by 



Am^g = Ami;, 



(3) 



as required by the CPT theorem [8] . Note that contribu- 
tions to the hamiltonian proportional to the unit matrix 
generate no oscillation effects, but they may nonetheless 
be relevant to stability and causality of the underlying 
theory [41]. 

The Lorentz-invariant hamiltonian ho can be diagonal- 
ized using a 6x6 unitary matrix U, 



{ho)A>B> = Ua'aUb'b{^o)ab, 



(4) 



AB 



where primed indices indicate the diagonal mass basis. 
The absence of mixing between neutrinos and antineu- 
trinos implies the mixing matrix is block diagonal, 



U = 



Ua'b 

U-,'-t 



with vanishing 3x3 ofF-diagonal blocks, 

Ua'b - U^'h = 0. 



(5) 



(6) 



Since the mass matrices for neutrinos and antineutrinos 
are related by complex conjugation, we also have 



Ua'b — U*ih- 



a'b' 



(7) 



The diagonal 3x3 blocks of ho can therefore be written 
as 



{ho)ab = {ho) 



ab / J 

a' = l,2.3 



Ul'JJa'bEa', 



(8) 



where Ea' are the usual three neutrino eigenenergies. In 
what follows, we assume that these three eigenenergies 
are nondegenerate. Note that this implies there are three 
twofold degeneracies in the full 6x6 hamiltonian (/io)ab- 
The Lorentz- violating term 5h'm Eq. (1) can be written 
in the form 



6h 



Shab 5h^i 
5hab Shab 



(9) 



For Lorentz- violating operators of renormalizable dimen- 
sion, the upper-left diagonal block takes the form 



Shab 



[(aL)"Pa-(ci)"Vp;3],, (10) 



and leads to mixing between neutrinos, where the 
neutrino energy-momentum 4-vector is denoted Pa = 
{E, -p) Ri E{1, -p) and (az,)^;,, {cl)"^ are complex co- 
efficients for Lorentz violation [13] . Hermiticity implies 



6hab = Shl^, 
which imposes the conditions 

{clTJ = {c.Tb^*- 



(11) 



(12) 



Similarly, the lower-right diagonal block of 5h produces 
mixing between antineutrinos, 

^Kb = ■^[{aRTPcc-{CR)°''^Pc,Pi3\-j^ 

= ^[-{aLrpo.-{cLr'Po.pp\i,- (13) 

The off-diagonal 3x3 blocks of 6h, Sh^j, and Shf,a, lead 
to neutrino-antineutrino mixing, an unconventional ef- 
fect. These blocks obey the hermiticity condition 



Sh^l = Shl^ 



(14) 



and can be written as [13] 



Sh 



ab 



iV2{e+)4r%-H-]^-^, 



6h-ab = zx/2(e+)*[rt/3-ff1,, 

= i^/2ie+)l[r% + HX-^. (15) 

In these equations, the complex coefficients for Lorentz 
violation g^^ and H"^ obey the relations 



rj 



^ba ifba ' 



^ab = -^ba = ^ba ' 



(16) 



The complex 4-vector (e+)a = (0, — e+) represents the 
helicity state. Introducing the local beam direction 
and other unit vectors associated with local spherical co- 
ordinates as 



Cr = (sin 9 cos 4>, sin 9 sin (f>, cos 6*) , 
ee = (cos 6 cos c^, cos 6 sin <p, — sin 9) , 
= (— sin(/>, COS0, 0), 

the 3-vector e+ can be expressed as 

1 



V2 



(ee + ie^) 



(17) 



(18) 



The coefficients (0^)"^ and H^j^ have dimensions of 
mass, and each taken alone leads to oscillation effects 
that are energy independent. In contrast, the coefficients 
{cL)'ab ^n.'i H'ai ^^"^ dimensionless, so their effects nat- 
urally scale with energy. Note, however, that combina- 
tions of coefficients can produce involved energy depen- 
dences, including mimicking conventional mass terms via 
the Lorentz- violating seesaw mechanism [13, 14, 17]. 
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The coefficients (ol)"}, and g"^ control CPT-odd ef- 
fects, wliilc {cL)2b ^^"^ ^ab govern CPT-even ones. Con- 
sequently, CPT symmetry holds when {aL)2b and g^^ 
vanish, and we find the oscillation probabilities obey the 
relationship 



(CPT invariance). (19) 



The CP symmetry may nonetheless be violated, so the 



relation 



= Pp_^^pj_ may fail. Further discussion of 



CP and CPT tests is provided in Sec. Ill C below. 

All the coefficients discussed here are taken to be 
spacetimc constants, so that translational symmetry and 
energy-momentum conservation hold. If the Lorentz vi- 
olation is spontaneous [3], which may be ubiquitous in 
effective field theories [42], then the coefficients can be 
understood as expectation values of operators in the fun- 
damental theory. Under these circumstances, requiring 
constancy of the coefficients is equivalent to disregarding 
soliton solutions and massive or Nambu-Goldstone (NG) 
modes [43]. When gravity is included, the NG modes can 
play the role of the graviton [44] , the photon in Einstein- 
Maxwell theory [45], or various new spin-dependent [46] 
or spin-independent [47] forces. The presence of gravity 
can also produce additional Lorentz-violating effects on 
neutrino oscillations [7, 48]. 

When neutrinos propagate in matter, the resulting for- 
ward scattering on electrons, protons, and neutrons can 
affect neutrino oscillations [49] . In the rest frame of the 
matter, this adds to the effective hamiltonian {hes)AB 
terms equivalent to CPT-odd coefficients given by [13] 

iaL,^«)l^. = {aL,^s)rr = -GFnn/V2, (20) 

where rig and n„ are the number densities of electrons and 
neutrons in the matter and Gf is the Fermi coupling con- 
stant. For example, in neutrino-oscillation experiments 
with long baselines, the propagation is over compara- 
tively long distances in the Earth's crust. In this case, the 
densities rig and n„ can be taken equal and constant to 
a good approximation, with v^Gi^rie — 2.1 x 10~^^ GeV 
~ (940 km)~^. In the perturbative analysis of Lorentz 
violation that follows, any matter effects can be taken 
as part of the unperturbed hamiltonian (/io)ab- In sit- 
uations where mass oscillations dominate, the matter ef- 
fects could alternatively be treated as perturbative and 
included as part of the Lorentz-violating term 6hAB- 



B. Perturbation series 

In this subsection, we use standard techniques of time- 
dependent perturbation theory to derive a perturbative 
series for the transition amplitudes. We treat the hamil- 
tonian component Sh describing Lorentz and CPT viola- 
tion as small compared to 1/L. 



The time-evolution operator S{t) is written in the form 

S{t) = e-'^'^«' 

= (e-*''<=«*e*''°*)S'("^(t) 

= ^("^(t) + S'^^\t) + S'^^Ht) + • • • , (21) 

where S^""^ is the nth-order perturbation in 6h. The con- 
ventional term is given by 



5(0) = 



ihot 



(22) 



The higher-order terms are obtained using the integral 
relation 



.+ [ dti{-i)Ah{ti) 
Jo 



t rt2 

dt2 / dti{-i)l\h{ti){-i)l^h{t2) + ■■■ , (23) 

^0 



where 



A/i(t) 



— ihot 



She 



ihcft 



(24) 



The integrals (23) can conveniently be performed in 
the mass-diagonal basis. To second order in Lorentz- 
violating coefficients, the results take the form 

= -^t5hA'B'r%,{t), 

= -ht'Y.^hA'C'5hc'B'T%c'it)- (25) 
c 

All sums over flavor indices are written explicitly 

throughout this work. The time dependence is contained 

in the factors r^?, r^,'^,, and t^a''b'C'- '^^^ zeroth-order 
factor is the usual expression 

The first-order term is given by 



(26) 



exp(— ii^g/t) / dii exp(— zA^/B/ti) 



= < 



exp(— ii^B'f), 

exp(— iii^^/t)— exp(— iS^/t) 



Ea' = Eb' 



otherwise, 



(27) 



where 



A^'B' = Ea' — Eb' 



(28) 

are the standard eigcncncrgy differences. The second- 
order factor is given by the integral 

^A'B'C'W = ^ex.p{-iEB>t) 

X / dt2 dtiexp{-iAA'C'ti)exp{-iAc'B't2) 
Jo Jo 

exp{-iEB't), Ea' = Eb' = Ec, 

(1) _ (1) (1) _ (1) 

I ^ ^^^Co't' = ^ ^^C^t' ' otherwise. 

(29) 
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We give two expressions in the last case so that expres- 
sions for the hmiting cases Ej^/ = Ec and E^' = Ec can 

readily be extracted. Note that both T^}g, and 7-^/3/(7/ 
are dimensionlcss functions of EA't that are totally sym- 
metric in mass-basis indices. 

Transforming to the flavor basis, the Lorentz-invariant 
transition amplitude is found to be 



c(0) 



J2ul'AUA'Be- 



(30) 



This leads to the usual oscillation probabilities for the 
Lorentz-invariant case of massive neutrinos. At first or- 
der, we choose to express the transition amplitude in the 
convenient form 



-it^{MAB)cD6hcD, 

CD 



(31) 



where the factors 

{M^Il)cD{t) = \j\HSfc{ti)stl{t-h) (32) 

= Y.''A'B'{mi'AUA'cU%,nUB'B 
A'B' 

depend on the energy and baseline of the experiment and 
also on the conventional masses and mixing angles. For 
given mass spectrum and mixing angles, these factors 
determine the sensitivity of an experiment. They are in- 
dependent of the direction of the neutrino propagation 
and of the coefficients for Lorentz violation. As a result, 
they remain unchanged as the Earth rotates. The quan- 
tity defined in Eq. (31) is a linear combination of 
these factors and the Lorcntz-violating perturbation 5h. 
It plays a key role in the expressions for the oscillation 
probabilities derived in the next subsection. Note that 
'H^AB reduces to SHab in the limit of negligible mass mix- 
ing. 

The second-order result for the transition amplitude 
can be written in a similar form. We define 



_1.2o,(2) 
2'' '^AB 



2*- 



^ {M%)cDEF5hcD5hEF,m 



CDEF 



where the experiment-dependent factors 

{M%)cDEF{t) 



ABICDEF 

Jo JQ 

(2) 



dt2 I dhS^^l{ti)S^^l;{t2-ti)SPs{t-t2) 

= E r'A,'B'C'it)Ul'AUA'cU^>DUc'EU*s,FUB'B 

(34) 



A'B'C 



again determine the combinations of coefficients relevant 
for oscillation effects. In this case, however, the quantity 



n'-Ali defined in Eq. (33) is a quadratic combination of 
coefficients. This leads to sidereal variations in neutrino 
oscillations at higher multiples of the Earth's rotation 
frequency. Note that H^ab reduces to {Sh'^)AB in the 
limit of negligible mass mixing. 



C. Probabilities 

Using the above results for the transition amplitudes, 
we can derive the oscillation probabilities. At zeroth or- 
der, the transition amplitudes are Lorentz invariant and 



take the usual block-diagonal form, S r 

^ ' ab 



(0) 



0(0) 
'-'ah 



0. So 



the zeroth-order probabilities for neutrino- antineutrino 
oscillations vanish, 



p(0) ^ p(0) ^ p 



(35) 



Since CPT is conserved whenever Lorentz symmetry 



holds [8], we have 5^°^ = S^^J . This implies 



(0) 



p(0) 



p(0) 



(0)1 2 

ab I ' 



(36) 



which leads to the usual results for Lorentz-invariant 
oscillation probabilities in terms of mass-squared differ- 
ences and mixing angles. 

The full mixing probability is given by 



?(0) 



, c(l) , q\^) I 



MS 



(2) 



AB 



(37) 



At second order in 5h, this gives 



p(0) 

p(l) 



p(2) 



■'as I 



I CV"; |2 



2Rc ((^("^ 



AB/ 



Ai)) 

'-ABJ 

(2)N , 

ab) ^ 



1^ 



(1) |2 



= -t'Re{{SZrn'IL)+t'\n%\'. (38) 

These equations involve the 6-dimcnsional space spanned 
by A. They can be decomposed into oscillation probabil- 
ities expressed in terms of the neutrino and antineutrino 
subspaces spanned by a and a. 

At first order, a short calculation shows that the prob- 
abilities can be written 



p(i) 

p(l) 
p(l) 



= 2tIm((^(°))*H« 
= 2tIm((4°))*H«: 
= = 0. 



(39) 



No neutrino-antineutrino mixing occurs because s'^^^ = 



S-^? = 0. Only the combinations Ti^V and H-V obtained 

ab •/ ab ab 

from the definition (31) contribute to these probabilities. 



(1) 
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Explicitly, we find 



to 



n 



ab 

(1) 



n 



(1) 

ab 



ab 



cd 

cd 
cd 



(40) 



Although nl^J and H]^^ are absent from the first-order 
probabilities, we include their expressions here because 
they enter the second-order probabilities below. Since 
the first-order results involve the diagonal blocks of 6h, 
only the coefficients {aL)ab a^nd (cl)"^ play a role. 

Decomposing the results (38) reveals that the second- 
order probabilities are 



cd 



p(2) 
p(2) 
p(2) 
p(2) 



= -t2Re((sW)*Hi- 



ab 



(41) 



H 



(2) ^ 
ab 



n 



(2_) 
ab 



^(2) _ 



n 



(2) 
ab 



cdef 

+ Yl (■^ab)cdef^Kd^f^ef, 

cdef 

Y:iM^Shd,f-sh,dSh,f- 

cdef 

+ Y.^M'S)s<ief-ShsdSKj, 

cdef 

J2(-^'kb)cdefShcdShef 

cdef 

cdef 
cdef 

+ T.iM'S)cd,f-SKdSh,f-. 

cdef 



(42) 



(2) (2) 

Note that Hl,u and H r are absent from the second- 

ab 

order probabilities but are included here for complete- 
ness. This implies that no cross terms between 5hab, 
Shg^f,, and 5h^i = 5h^^ appear at second order. In par- 
ticular, all appearances of the coefficients g^^ and H^i 

arise as squares or as quadratic products with each other, 
without accompanying factors of (aL)a6 i'^i^)'ab ■ 



all 



where S 



(0) 



S. 



(0) 



is used to simplify the last 



two. The probabilities Pi 



(2) 



and P^^ 



(2) 



include both 



leading-order contributions from H.)^^ and H.-^ as well 
as higher-order contributions from the combinations H^^^ 
and 70. Also, nonzero mixing between neutrinos and 

ab ^ ..... 

antincutrinos appears, giving sensitivity to the linear 
combinations T-S^') and Ti}^) . This shows that the domi- 

ab ab ^ 

nant effects of the coefficients 'g^ and appear only 

at second order. Moreover, (ol)"}, or (cL)"f play no role 
in neutrino-antineutrino mixing at this order. 



Explicit expressions for W^^^ and W^^^ can be obtained 
by decomposing the quadratic combinations defined 

in Eq. (33). The structure of the factors {M\q)cdef 
and the form of the mixing matrix JJ reduce the num- 
ber of terms that contribute. In particular, we find 
that {M.\q)cdef vanishes unless the index pairs {AC}, 
and {DE} lie in the same subspace. This leads 



(2) 



III. COEFFICIENTS {aLTai, AND {cl)Z 



In this section, we consider effects originating from 
the coefficients (0^)"^ and [chf^ ■ These contribute 
only to neutrino-neutrino mixing and to antineutrino- 
antineutrino mixing. We focus here on the dominant 
signals, which arise from the first-order oscillation prob- 



abilities P, 



(1) 



and P'^J_,jj^ given in Eq. (39). 



The theoretical analysis presented in the previous sec- 
tion applies to any scenario involving neutrino propaga- 
tion. However, the key experimental signals are different 
for beam experiments, solar-neutrino studies, and cos- 
mological observations. For definiteness in this work, we 
provide rcsTilts in the context of beam experiments. 

Section III A establishes the key expressions describing 
sidereal variations in the oscillation probabilities. Illus- 
trative examples involving all three generations arc pro- 
vided in Sec. IIIB, while the two-flavor case is considered 
in Sec. Ill B 3. The construction of CP and CPT asym- 
metries to characterize the effects is considered in Sec. 

inc. 
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A. Sidereal variations 

The combinations of coefficients for Lorentz and CPT 
violation entering the nonzero first-order probabilities 
(39) are controlled by the experimental factors {■M''^^)cd 
and (A^g^^)gj entering Eq. (40). These factors can be cal- 
culated from Eq. (32). The time t can be set equal to the 
baseline distance L because any difference between the 
two is small and leads only to suppressed higher-order 
corrections. For given values of E and L, the nine com- 
plex constants (A4l^ij)cd determine the coefficient com- 
binations relevant for v ^ v mixing, while the nine 
complex constants {A4^}^)gfi determine those for 9 i-^ 9 
mixing. If CP is conserved in the usual mass matrix, 
as occurs in the two-generation limit, then the mixing 
matrices are real and obey Uab = t^aS- We then find 
{Ml^fj)cd = (A^g^^)gj, so a single set of nine constants 
determines the experimentally relevant combinations for 
both neutrinos and antineutrinos. 

The specific combinations of (ol)";, and (cl)"^ appear- 
ing in the transition probabilities can be found by consid- 
ering the forms of the relevant blocks of the hamiltonian. 
Eqs. (10) and (13). In this context, it is convenient to 
define the linear combinations 

(SL)«6 = ^(Ml;^)cd(aL)?d, 

cd 
cd 
cd 

= E(^S)c-d-M?|. (43) 

cd 

Ttic^sc; arc experiment-dependent combinations of the 
fmidamcntal coefficients {aL)ab ^^d {cL)ab ■ Using 
(aL)afe: (afl)a6, (cl)^^ , and {cr)^^ , the combinations 
W^^^ and Wgj^ controlling the first-order probabilities can 
be written in a form that mimics the hamiltonian pertur- 
bations dhab and (J/igb, 

^is = |[(«K)"Pa-(?«)"Vw],5. (44) 

This form reveals the explicit 4-momentum dependence 
of the transition probabilities. 

The momentum dependence implies that the mixing 
behavior can depend on the direction of neutrino prop- 
agation. For Earth-based experiments, the source and 
detector rotate at the sidereal frequency uj(^ ~ 27r/(23 h 
56 min), which can induce sidereal variations in the os- 
cillation probabilities. Since the first-order probabilities 
are linear in Ti.'^V or 7i*:V, each of which has both linear 
and quadratic terms in the 3-momentum, sidereal varia- 



tions controlled by the coefficients (a^)";, and (ci,)"f can 
occur at the frequencies and 2loq. 

To display explicitly these variations, a choice of in- 

ertial frame must be specified. By convention and 
convenience, the standard incrtial frame is taken as a 
Sun-centered celestial-equatorial frame with coordinates 
{T,X,Y,Z) [9, 50]. The Z axis of this frame is di- 
rected north and parallel to the rotational axis of the 
Earth. The X axis points from the Sun towards the ver- 
nal equinox, while the Y axis completes a right-handed 
system. The origin of the time coordinate is chosen as 
the vernal equinox 2000. The Earth's rotation causes the 
neutrino 3-momentum to vary in local sidereal time 
at the frequency cjq in the Sun-centered frame, unless it 
happens to lie along the Z axis. 

For neutrino- neutrino mixing, we can display explicitly 
the sidereal variation by expanding H^^^ as 



+ {Ai^'>)abSmiJ^T(s + (Ai^'^) ab cos CJqT(s 

-h(e(i))„bsin2wel^e + {Bi'^)ab cos 2us,T^. (45) 



Suppose the neutrinos of interest are emitted in a defi- 
nite direction relative to the Earth, perhaps as a neutrino 
beam from an accelerator. Let the vector {N^ , , N^) 
represent the propagation direction in the Sun-ccntcrcd 
frame at local sidereal time = 0. We can write this 
vector in terms of local spherical coordinates at the de- 
tector. Denote by x the colatitudc of the detector. In- 
troduce at the detector the angle between the beam 
direction and vertical, and also the angle </> between the 
beam and east of south. The components of the vector 
can then be written as 



TV-^ = cos X sin 6 cos + sin x cos 6, 

= sin 6 sin (f), 

= — sinxsin0cos0 -|- COSXCOS0. (46) 



Using these expressions, the amplitudes in the expansion 
(45) are specified in terms of coefficients for Lorentz vio- 
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lation in the Sun-centered frame as 

+ \{l-ZN^N^)E{cLfa^, 

+2N''E{Zl)1^ + '^N''E{Zl)1^ 
-2N''N^E{Zi.)l^ - 2N^N^E{cl):,^, 



-{N^N^-N^N^)E{cl)1'', 
iBi')U = -2N^N^Eic,)^,y 



l{N^N^-N^Ny)E{{cj^)^,^ 



(47) 



The form of the above expansion matches that used in 
Ref. [16] in the context of short-basehne neutrino exper- 
iments. 

The sidereal variation in H^^^ leads to a correspond- 
ing variation in the probabilities. We parametrize these 
variations as 



p(i) 



2L 



{P^^'>)ab + (Pi^)„6Sinc^eTffi + {P^^!)abCOSUJ(sTe 



+(p£^)„6sin2a;eTe + {P^'JUcos2u<sTq, 

(48) 



where 



{P^U = Im((S(°))*(C«)„,), 

(P«)„, = Im((5f,V(4^))„,), 

(P«)„b = Im((Si°))*(4i))„b), 

{P^Wt, = Im((5(°))*(B(i))„6), 

{P^%b = lm{{S^^^nBP)ab) (49) 

are combination of coefficients for Lorentz violations. 

Note that the sidereal amplitudes {P^^^)ab, {PA^)ab, 



(■Pa?)"'" iPB,')ab, and {PB^)ab are tiny, with size deter- 
mined by {aL)ab and (cl),^ • The expression (48) reveals 
that the experimental sensitivity to perturbative Lorentz 
and CPT violation increases with the baseline L. 

For antineutrino-antineutrino oscillations, analogous 
results hold. We can expand in the form (45), 
replacing the indices {ab} with {56} . The amplitudes 
again take the form (47), but with the_ substitutions 
(Si) (clr), (cl) (cr), {ab} {ab}. Similarly, 
the sidereal variation in the probabilities can be written 
as Eq. (48) by replacement of the indices. 



7(1)^ 



B. Illustrations 

In this subsection, the first-order perturbative formal- 
ism derived above is applied to several illustrative situa- 
tions. Following some comments about the methodology, 
we consider an example involving mixing of all three fla- 
vors of neutrinos and then discuss the limiting case of 
two flavors. 



1. Methodology 

Since the effects from Lorentz and CPT violation are 
perturbative, an explicit analysis must specify the con- 
ventional mass spectrum and mixing angles. In the stan- 
dard three-neutrino massive model [2] , the usual 3 x 3 ef- 
fective hamiltonian {ho)ab for neutrino-neutrino vacuum 
mixing appears as the upper-left block of Eq. (2). It can 
be written as 

1 



{ho)ab = ^^^Ib 



^J^^a'aUb'b^mlt,,, 

a'b' 



(50) 



where Am^/,/ is the diagonal mass matrix. Only 
two mass-squared differences contribute to oscillations. 
Without loss of generality, we can therefore express the 
diagonal mass matrix as 



m 



a'b' 



'Q 

Am| 
,0 At- 



(51) 



For a normal mass hierarchy, the quantity Atoq is the 
smaller of the two mass-squared differences. It is of 
particular relevance in situations involving low-energy 
oscillations such as solar- neutrino measurements. The 
larger mass-squared difference Am\^^ plays a central role 
where mixing of high-energy neutrinos occurs, such as 
atmospheric-neutrino experiments. The mixing matrix 
Ua'b can be written as 





( Cl2 


-Sl2 




Ua'b = 1 


.Sl2 


Cl2 


( 




I 









(52) 



where Cij = cos( 



Sij = sinOij, and d is the CP- 



violating phase. For antineutrino-antineutrino mixing, 

the conventional effective hamiltonian for vacuum mixing 
is obtained by complex conjugation of the above results, 
which is equivalent to changing the sign of 5. 

In the presence of matter, the 3x3 effective hamil- 
tonian {ho)ab for neutrino-neutrino mixing acquires an 
additional term from Eq. (20) and becomes 

1 

2E' 



{ho)ab 



bei 



(53) 
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where the index e refers to the electron-neutrino flavor. 
The additional term changes the eigenvalues of (/io)a6 ^md 
hence the explicit values of the components of the overall 
mixing matrix. Note that the expression (53) is strictly 
valid only in the rest frame of the matter. For the specific 
scenarios considered below, this frame comoves with the 
Earth as it rotates on its axis and revolves about the 
Sun. These motions are nonrelativistic, however, so the 
comoving frame can be identified with the Sun-centered 
frame of Sec. Ill A to an accuracy of parts in 10**. We 
adopt this identification in what follows. 

In performing an analysis for Lorentz and CPT viola- 
tion, the appropriate methodology depends on the loca- 
tion of the experiment in L-E space and on the presently 
unknown value of ^13. Consider, for example, three hy- 
pothetical beam experiments with the same long base- 
line L ~ 200 km but with differing energies i?i ~ 10 
MeV, E2 ~ 1 GeV, £3 ~ 100 GeV. The first lies in a 
region where mass oscillations involve all three flavors, 
so the treatment of Lorentz and CPT violation requires 
the three-flavor formalism of the previous section. The 
same is true of the second experiment for large 6*13. How- 
ever, if f?i3 is small or zero, then significant mass mixings 
in this second experiment involve only two generations. 
Lorentz- violating effects within these two generations can 
then be studied using a two-flavor limit of the previous 
section, while effects involving the third flavor are well 
described using the procedures for negligible mass mix- 
ing presented in Ref. [16]. For the third experiment, no 
significant mass mixings occur and so the methods of Ref. 
[16] are applicable for Lorentz and CPT violation in all 
flavors of oscillations. 

High sensitivity to operators for Lorentz and CPT vi- 
olation of mass dimension three, such as those controlled 
by the coeflicients (aL)ab, can be achieved in experiments 
with long baselines L. Similarly, high sensitivity to op- 
erators of mass dimension four, such as those governed 
by {cl)"^ , can be obtained via long baselines L, high en- 
ergies E, or both. At present, most existing or planned 
long-baseline experiments have baselines L ^ 200-1500 
km and energies ~ 1-10 GeV and hence lie in a region 
of L-E space analogous to that of the second hypothetical 
experiment above. 

For illustrative purposes, we consider here a variety of 
beam experiments involving long baselines L and seeking 
Vf, appearance in z/^ beams or studying disappearance. 
Existing experiments in this category include KEK to 
Super-Kamiokande (K2K) with baseline L ~ 250 km [51], 
the MINOS far detector with baseline L ~ 750 km [52], 
and the Oscillation Project with Emulsion- Tracking Ap- 
paratus (OPERA) [53]. The latter has essentially iden- 
tical baseline L ~ 750 km to the Imaging Cosmic and 
Rare Underground Signals experiment (ICARUS) [54]. 
The Fermilab E929 experiment (NOz/A) [55] with base- 
line L ~ 800 km is currently under construction, while 
the Tokai to Kamioka (T2K) experiment [56] with base- 
line L ~ 300 km has recently begun data taking. Other 
experiments with even longer baselines are under con- 



sideration, including one at the Deep Underground Sci- 
ence and Engineering Lab (DUSEL) [57] using a neutrino 
beam from Fermilab and baseline L ^ 1300 km, and the 
Tokai to Kamioka and Korea (T2KK) experiment [58] 
using the same neutrino beam as T2K but with baseline 
L ~ 1000 km. All of these experiments have excellent 
sensitivity to perturbative Lorentz and CPT violation. 

For definiteness in what follows, we consider two ex- 
plicit scenarios. In the first, discussed in Sec. IIIB2, we 
take a comparatively large value of ^13 and consider the 
effects of Lorentz and CPT violation for studies of ap- 
pearance. For this situation, mass mixing involving all 
three fiavors occurs and so the full formalism of the previ- 
ous section is appropriate for the analysis. In the second 
scenario, considered in Sec. TUB 3, we suppose 6*13 is neg- 
ligible and investigate the effects of Lorentz and CPT vi- 
olation on i/p disappearance. For this case, only ^ 
involves significant mass mixing and so a two-flavor limit 
of the previous section can be applied. 

We can also identify two interesting regions of L-E 
space that could benefit from the development of new 
experiments. The first is the region of long baselines 
L > 200 km with low energies E < 100 MeV. This is of 
particular interest if ^13 ~ 0°, since it provides an op- 
portunity for clean studies of three-flavor mixings that 
are otherwise challenging to perform. A comparatively 
intense source is needed due to the long baseline and 
the cross-section falloff with energy. One possibility is 
a setup similar to the Kamioka Liquid Scintillator An- 
tineutrino Detector (KamLAND) [59] but with a single 
source, directional sensitivity, or both. For longer base- 
lines, a beta beam [60] may be an interesting option. A 
low-energy beta beam has been studied in the context 
of short-baseline experiments [61]. Another possibility 
might be an intense pulsed neutrino beam such as that 
proposed for a neutrino facility at the Spallation Neutron 
Source (z^-SNS) [62] and lying in the 10-50 MeV energy 
range. 

The second interesting region lies at high energies 
E > 100 GeV, even for comparatively short baselines 
L < 10 km. Neutrino beams at these energies have 
been used in the Neutrinos at the Tevatron (NuTeV) [63] 
and Chicago-Columbia-Fermilab-Rochester (CCFR) [64] 
experiments. For studies of the coefficients (cl)"^ for 
Lorentz and CPT violation, a high energy compensates 
for a shorter baseline and so oscillation experiments in 
this region could have a competitive reach. Since mass 
mixing is negligible, the methods of Ref. [16] are appli- 
cable for analyzing Lorentz and CPT violation in this 
case. 



2. Example: Ve appearance 

Consider first searches for Lorentz and CPT violation 
via Ue appearance in a beam, within the assump- 
tion of a comparatively large ^13. This requires a full 
three-flavor analysis. We incorporate effects from matter- 
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Experiment 


K2K 


MINOS 


OPERA, 
ICARUS 


NOz/A 


T2K 


DUSEL 


T2KK 




-0.05 


-0.10 


-0.01 


-0.17 


-0.16 


-0.08 


-0.11 




0.84 


0.63 


0.88 


0.38 


0.48 


0.38 


0.41 




-0.13 


-0.24 


-0.02 


-0.46 


-0.46 


-0.35 


-0.39 




0.00 


0.00 


0.00 


0.00 


0.00 


0.01 


0.01 


RoiMi'J)^^ 


-0.05 


-0.09 


-0.01 


-0.15 


-0.15 


-0.09 


-0.11 


n / A A{i-)\ 

Re{MliJ)i^T 


-0.01 


-0.02 


0.00 


-0.01 


-0.01 


0.04 


0.03 


Re (Mia) re 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


Ro(Mi]})-r„ 
' •j*-' \^>'*^ie/j. / T fj, 


-0.05 


-0.08 


-0.01 


-0.14 


-0.14 


-0.06 


-0.08 


Re(M9-]]TT 


-0.01 


-0.01 


0.00 


0.00 


0.00 


0.03 


0.02 




-0.08 


-0.06 


-0.02 


-0.03 


-0.08 


0.10 


0.07 


lm{Mi]I)ei. 


-0.39 


-0.62 


-0.44 


-0.61 


-0.43 


-0.34 


-0.32 




-0.24 


-0.24 


-0.06 


-0.21 


-0.27 


0.17 


0.12 




0.00 


0.01 


0.00 


0.02 


0.03 


-0.04 


-0.03 




-0.07 


-0.06 


-0.02 


-0.02 


-0.02 


-0.05 


-0.05 


Im (A4efi jpT 


U.Ui 


n no 
0.U2 


U.OO 


U.UI 


n no 
U.08 


n no 
-U.Oo 


-U.U6 




0.00 


0.01 


0.00 


0.02 


0.02 


-0.05 


-0.04 




-0.06 


-0.05 


-0.02 


0.00 


0.00 


-0.02 


-0.02 


Im {M^]})^r 


0.01 


0.02 


0.00 


0.06 


0.07 


-0.12 


-0.10 


numerical values of experiment-dependent dimensionless factors (A^i)2)cd 



MINOS, OPERA, ICARUS, NOi^A, T2K, DUSEL, and T2KK. Numerical values are computed via Eq. (32) adopting the 
parameters (54), using estimated beam baselines L and neutrino energies E for each experiment. Within this approximation, 
the antineutrino factors (A^^ji* )cd are identical. 



induced mixing via Eq. (53), and we adopt explicit pa- 
rameter values for vacuum mixing compatible with the 
observed three-flavor oscillations in solar, atmospheric, 
reactor, and accelerator experiments [2], 

Am| ~ 8.0 X 10-^eV^ 

Aml^ 0, 2.5X 10-''cV^ 
012 ^ 34°, 923 ^ 45°, 6»i3 12°, 6 :^ 0°. (54) 

Together with the matter effects, the values (54) deter- 
mine the linear combinations of coefficients {aL)2b ^^'^ 
(^i')af relevant for any given experiment. Note that 
both magnitudes and signs arc significant. For exam- 
ple, the above values hold for a normal mass hierarchy. 
Our general expressions for oscillation probabilities are 
valid for any magnitudes and signs, but the numerical 
results for the illustrative examples below assume the 
specific choices (54) . For a comprehensive exploration of 
the space of coefRcients for Lorentz and CPT violation, 
distinct analyses of the same data must be performed 
for each acceptable choice of parameter values and must 
be reported as such. Note also that the approximation 
^ ~ 0° made in Eq. (54) implies that there is little or 
no CP violation in standard oscillations, although per- 
turbative CP violation from Lorentz and CPT violation 
can still appear. 

The anisotropics introduced by nonzero coefficients for 
Lorentz violation can lead to sidereal variations, which 



are the primary signals of interest here. The sidereal de- 
composition of the probability takes the form (48) with 

{ab} = {e^}. The four amplitudes {P'"2l)efj,, (-P4j)eM, 

(Pg^^)e^, and {Pg^)ep, are linear combinations of coefR- 
cients for Lorentz violation given by Eqs. (47) and (49), 
and they can be measured by studying the variations of 
the neutrino mixing with sidereal time. The effects from 
the combination {P^^^)cfj, are more challenging to detect 
experimentally because they have no accompanying time 
variation. 

The amplitudes of the sidereal-variation probabili- 
ties depend on the quantities (ol)"/* and (cl)"^, which 
are linear combinations of fundamental coefRcients for 
Lorentz and CPT violation, 

cd 

(<^L)t = T.(^eJU{cLrJ. (55) 
cd 

The results of an experimental analysis can therefore be 

expressed in terms of the fundamental coefhcients (aL)cd 
and (cl)"^ by calculating the relevant complex factors 

{■M.e]})cd for the given experiment. 
Table I presents approximate numerical values of 

the real and imaginary parts of these factors for the 
eight long-baseline experiments K2K, MINOS, OPERA, 
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Experiment 



K2K 



MINOS 



OPERA, ICARUS 



Amplitude 



p(i) 



p(l) 



p(l) 



p(l) p 



(1) p{l) p{l) 



Re (Si)^ 

Re(?L)f/ 
Re(?L)f/ 

YZ 



p(l) 


p(i) 
Oi, 


p(i) 


p(l) 


0.0 


0.0 






0.0 


0.0 






-1.0 


0.6 






0.6 


1.0 







Re (ci 
Re(cz,), 

Im (az,)e 
Im (ol), 
Im {cl)1 

Im (cL)e 



) 

Im (ol)^^ 
Im (Hi)f/ 

\TY 

Im (cz,)f/ 
Im (Hi)f/ 
Im (cL)f/ 
Im (cz,)^/ 
Im (cL)r/ 
Sensitivity 



-0.1 0.0 

0.0 0.1 

0.4 0.0 

0.0 -0.4 - 

0.0 0.1 
0.0 -0.1 

0.0 0.0 

0.2 0.0 

0.0 0.0 

0.1 0.0 
0.0 -0.1 - 
-0.3 0.0 
0.0 0.0 

0.0 -0.1 
0.0 0.1 

0.0 0.0 

- -0.1 0.0 
0.0 0.0 

< 8 X 10"^^ 



0.0 0.1 
0.1 0.0 
0.2 -0.3 - 
-0.3 -0.2 - 

0.1 0.0 

- -0.1 0.0 
-0.1 0.2 

0.0 -0.1 

0.2 0.1 

0.1 -0.1 - 

-0.1 -0.1 - 
-0.4 0.5 
0.2 -0.3 - 

- -0.1 0.0 

0.1 0.0 
0.2 -0.3 - 

0.1 0.2 
-0.3 -0.2 - 
< 3 X 10"^^ 



-0.4 0.2 



0.2 0.4 



0.2 0.1 
-0.2 -0.1 

0.3 -0.4 



0.0 0.0 
0.0 0.0 
0.5 -0.3 - 
0.2 -0.1 - 

- -0.1 -0.1 

0.1 0.1 
0.2 -0.1 - 

- -0.1 0.3 
-0.1 -0.2 - 

< 3 X 10"^^ 



Experiment 



NOi'A 



T2K 



DUSEL 



T2KK 



Amplitude 



p(i) p(i) piJ-; p' 
^As ^Aa ^B, 



(1) 



p(l) p(l) p(l) p(l) 

-f^A, ^Ao ^Bs ^B^ 



p(l) pKi-l p 
^As ^Ac ^B 



(1) 



(1) 



Be 



Re(oi)^ 

Re(aL)r^ 
Re(?L)i;f 

\TY 
\XX 

\XZ 
Je/j. 
\YY 

Re (Cz,)ej^ 

Im(oi)f^ 
Im (Sl)!^ 

Im (cl)^J 

Im (ci,)f/ 
Im (cL)f/ 
Im(?£)f/ 
Im (cL)rM^ 



p(l) 

^A, 


p(l) 


-0.1 


0.0 


0.0 


0.1 


0.1 


0.0 


0.0 


-0.1 



py-) py- 

^B. ^Ba 



Re(£i) 
Re(ci)^ 

Re (CL)e 
Re (CL)e 



Sensitivity 



0.0 0.0 
0.0 0.0 
0.0 0.0 
0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.2 -0.2 - 
-0.2 -0.2 
-0.6 0.8 
0.4 -0.5 - 

- -0.2 0.1 
0.2 -0.1 
0.4 -0.5 - 

0.1 0.4 
-0.5 -0.4 - 



0.0 0.0 
0.0 0.0 
0.0 0.0 
0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.3 0.0 

0.0 -0.3 - 
-0.4 0.0 
0.0 0.0 

0.0 -0.1 
0.0 0.1 

0.0 0.0 

- -0.2 0.0 
0.0 0.0 



-0.1 0.0 

0.0 0.1 

0.1 0.0 

0.0 -0.1 - 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.1 0.0 

0.0 0.0 

0.1 0.0 

0.0 -0.1 - 
-0.3 0.0 
0.0 0.0 

0.0 -0.1 
0.0 0.1 

0.0 0.0 

- -0.1 0.0 
0.0 0.0 



0.0 0.0 



0.0 0.0 



0.0 0.0 

0.0 0.0 

0.1 0.0 



0.2 0.0 

0.0 -0.2 - 
-0.3 0.0 
0.0 0.0 

0.0 -0.1 
0.0 0.1 

0.0 0.0 

- -0.2 0.0 
0.0 0.0 



< 2 X 10 



-23 



< 7 X 10 



-23 



< 2 X 10- 



-23 



< 2 X 10 



-23 



TABLE II: Estimated amplitudes of sidereal-vaxiation probabilities {P^'')eij., (-PA^je/^, {PB')eij,, and {P^^')eix for appearance 



and (P«) 

experiments with p^- The numerical value is listed for the estimated contribution to each amplitude from the real and 

imaginary parts of the combinations (aL)"^ and (cl)"^ of fundamental coefficients for Lorentz and CPT violation. Values arc 
given to one decimal place, in dimensionless units for (5^)"^ and in units of GeV for (ci,)"^. A value of 0.0 indicates rounding 
to zero at this precision, while a dash implies the value is identically zero. The last row lists the approximate sensitivity in 
GeV of each experiment. 
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ICARUS, NOi^A, T2K, DUSEL, and T2KK. The en- 
tries are obtained assuming the parameter values (54) 
and incorporating effects of matter-induced mixing. The 
factors vary with the experiment, reflecting their depen- 
dence on the basehne L and the neutrino energy E. Since 
the experiments are clustered in a single region of L-E 
space, the corresponding numerical values for each factor 
are roughly comparable and so the four sets of (ol)"^, 
(ci)g^ obtained from Eq. (55) are roughly compara- 
ble combinations of the fundamental coefficients (ol)"},, 
(cl)^^. The table reveals that each experiment can mea- 
sure particular combinations of (fti)"^ and (cl)"^, which 
can then be used to constrain the coefficient space of 
{aL)2b ^^'^ i'^L)ab- The large number of coefficients for 
Lorentz and CPT violation and the limited number of 
observables for a given experiment imply that multiple 
experiments are needed to explore the entire coefficient 
space. In performing the analysis, it is of practical value 
to obtain an estimated maximal sensitivity to each indi- 
vidual component of the fundamental coefficients (oz,)";, 
and (ci)"^ in turn, by allowing only that component to 
be nonzero and using the data to constrain it. 

We can use Eqs. (47) and (49) to calculate estimated 
first-order sensitivities to (a^)"^ and (c^)"^ for each of 
these experiments. Table II lists the four amplitudes 

(P_J^)e^, (P^l)e^, (PB])e^, and (p£^)eM as explicit lin- 
ear combinations of the real and imaginary parts of 
(cLLjen ^"^^ The linear combinations typically dif- 

fer for each experiment and for each amplitude because 
they depend on the beam energy E and direction p. Note 
that the two experiments OPERA and ICARUS have ap- 
proximately the same baseline, orientation, and energy, 
so they can be listed together for our purposes. As an ex- 
ample, the table reveals that in GeV units the amplitude 
(P^^^)e^ for the K2K experiment is approximately 

(^il?)e/. ^ -0.1Re(Si)f^ + 0.4Re(?Oj/ 

-f0.1Im(5i)f^-0.3Im(?i)J^^. (56) 

Using Table I, all these combinations can be written in 
terms of the fundamental coefficients (ol)"}, and (cl)"^ 
for Lorentz and CPT violation. 

To obtain a crude estimate of the sensitivities for each 
experiment, we suppose that a 10% sidereal variation in 
the oscillation probability can be detected. This leads 
to a sensitivity of order 10%/2_L. The last row of Table 
II lists these values for each experiment. In conjunction 
with the other entries in Table II and with the factors 
listed in Table I, these values can be used to obtain the 
estimated first-order reach for any desired coefficient for 
Lorentz and CPT violation. 

Table II shows that long-baseline experiments have the 
potential to achieve extreme sensitivities to Lorentz and 
CPT violation in the v^. appearance mode if there 

is appreciable mass mixing arising from a comparatively 
large value of (?i3. Since some of the predicted effects 
are small, second-order effects may also play a role and 



may need to be incorporated in a comprehensive analysis 
of real data. If instead ^^13 is tiny or zero, studies of 
Lorentz and CPT violation in the — ^ appearance 
mode can be performed using the methodology of Ref. 
[16], as discussed in the previous subsection. 

In the event that the listed experiments run in an- 
tineutrino mode, the attainable reach can be estimated 
similarly. In the vacuum, the factors (mS^^-^^ are un- 
affected because the parameter values (54) imply CP in- 
variance, so the corresponding amplitudes of the sidereal- 
variation probabilities can be found by replacing (ol)"^ 

and (cl)"^ with (a^)"^ and (cr)"^. However, the con- 
tribution from mass-induced mixing in Eq. (53) enters 
with opposite sign, so the estimated amplitudes in Table 
II acquire corresponding changes. 

3. Example: disappearance 

Beams of u^^ also provide opportunities to search for 
Vfj^ disappearance. The probability of oscillating into 
other neutrinos is given by 

The correction introduced by {a,L)Zb ^'^^ (^^)af coeffi- 
cients is therefore given by 

Pl%., = -Pl%.^ = -2ilm ((5W)* W«), (58) 

where H^ill is defined in Eq. (44). Again, sidereal vari- 
ations can arise from the anisotropics introduced by 
Lorentz violation, and 7i|/J can be expanded in sidereal 
time according to Eq. (45). 

For v^^ disappearance experiments, the sidereal decom- 
position of the probability is given by Eq. (48) with 

{ab} = {wi}. The four amplitudes (Pjt^j)^^, {PaI)ij-iJ'^ 
of sidereal-variation probabilities 
are provided in Eqs. (47) and (49). Their exact expres- 
sions in terms of the combinations (Sl)^^ and (cl)^^ of 
fundamental coefficients for Lorentz and CPT violation 
take the form 

cd 

i^L)f, = Y^iM^^^UicLrJ- (59) 

cd 

If needed, the experiment-dependent complex factors 
(A^A*2)cd controlling these linear combinations can be ob- 
tained using Eq. (32). 

In scenarios with oscillations occurring primarily be- 
tween i/fj^ and i/r, the probability for i^^ disappearance 
can be well approximated by restricting attention to two- 
flavor vacuum mixing. This limit involves only one mass- 
squared difference and one mixing angle, and it offers an- 
other useful illustration of the general analysis given in 
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Experiment 


K2K 


MINOS 


OPERA, ICARUS 


Amplitude 


p(i) 


p(i) 


p(i) 


p(i) 


p(i) 

^As 


p(i) 

^Aa 


p(i) 


p(i) 

^Bc 


p(i) 

^As 


p(i) 

^Aa 


p(i) 

^B, 




Re(az,);Jv 


-0.5 


0.0 


— 


— 


-0.2 


0.2 


— 


— 


0.1 


-0.1 


— 


— 


Re(aL)^^ 


0.0 


0.5 


— 


— 


0.2 


0.2 


— 


— 


-0.1 


-0.1 


— 


— 




1.2 


-0.1 




— 


0.7 


-1.0 


— 


— 


-3.5 


2.0 


— 


— 




-0.1 


-1.2 




— 


-1.0 


-0.7 


— 


— 


2.0 


3.5 


— 


— 








Ii.il 


0.3 






0.2 


-0.1 






0.8 


0.5 


Re (cl)^^ 






0.0 


-0.3 






-0.2 


0.1 






-0.8 


-0.5 


Re(ci);f/ 


-0.1 


0.0 






-0.4 


0.6 






-1.4 


0.8 






Re(ci)^^ 






0.6 


-0.1 






-0.1 


-0.4 






0.9 


-1.6 


Re(ci)^T? 


0.0 


0.1 






0.6 


0.4 






0.8 


1.4 






Sensitivity 


< 8 X 10"^^ 


< 3 X 10"^^ 


< 3 X 10"^^ 



Experiment 


NOj/A 


T2K 


DUSEL 


T2KK 


Amplitude 


p(i) 

^As 


p(i) 

^Aa 


p(i) 


p(i) 


^As 


p(i) 

^Aa 


p(i) 

Bs 


p(l) 


p(i) 


p(i) 

^Aa 


p(i) 

Bs 


p(i) 


p(i) 


p(i) 

^Aa 


p(i) 


p(i) 

Be 




-0.1 


0.2 






-0.2 


0.0 






-0.5 


0.0 






-0.5 


0.0 






Re (aL)^T 


0.2 


0.1 






0.0 


0.2 






0.0 


0.5 






0.0 


0.5 








0.6 


-0.8 






0.2 


0.0 






1.0 


0.0 






0.8 


0.1 








-0.8 


-0.6 






0.0 


-0.2 






0.0 


-1.0 






0.1 


-0.8 






ReicL)^^"" 






0.2 


-0.1 






0.0 


0.1 






0.0 


0.2 






0.0 


0.2 


Re (cl)^^ 






-0.2 


0.1 






0.0 


-0.1 






0.0 


-0.2 






0.0 


-0.2 


Re{cL)^f 


-0.4 


0.5 






0.0 


0.0 






-0.2 


0.0 






0.1 


0.0 






ReicL)^^ 






-0.1 


-0.3 






0.1 


0.0 






0.5 


0.0 






0.4 


0.1 


Re{CL)^r 


0.5 


0.4 






0.0 


0.0 






0.0 


0.2 






0.0 


-0.1 






Sensitivity 


< 2 X 10"^^ 


< 7 X 10"^^ 


< 2 X 10"^^ 


< 2 X 10"^^ 



TABLE III: Estimated amplitudes of sidereal- variation probabilities {P^^)i^t, {Pac')!^'^' (^b1'')mti and {P^^)^t within the two- 
generation approximation. The numerical value is listed for the estimated contribution to each amplitude from the real parts 
of the fundamental coefficients (oi)™,- and {cl)"t for Lorentz and OPT violation. Values are given to one decimal place, in 
dimensionless units for (0^)^^ and in units of GoV for {cl)"t ■ A value of 0.0 indicates rounding to zero at this precision, 
while a dash implies the value is identically zero. The last row lists the approximate sensitivity in GeV of each experiment. 
Analogous results for possible experiments with antineutrinos can be obtained by replacing (ol)J^^ and (cz,)JJ^ with {anj^f and 



Sec. Ill A. In the event of tiny or zero ^13, it is also the 
relevant limit for the eight experiments considered above. 
We adopt this case as our second illustrative example. 

In the Lorentz- invariant two-flavor limit, an overall di- 
agonal term can be removed from the hamiltonian (/io)o6 
because it is irrelevant for oscillations. This gives a 2x2 
mass matrix of the form 



{ho)ab 






Ato|2 



u. 



(60) 



The flavor indices are now restricted to two generations, 
a,b,... = fi,T. The mixing matrix U depends on the 
mixing angle 623 according to 



Ua'a 



C23 -S23 
,S23 C23 , 



(61) 



Note that CP violation due to mass mixing is strictly 
unobservable in this limit. The single mass difference is 



then given by 



(62) 



For the explicit estimations in this subsection, we choose 
for Am^tjjj and 6*23 the values 



2.5 X 10-3 eV^ 

45°, 



(63) 



which are consistent with the three-flavor parameter val- 
ues (54). 

In the presence of Lorentz and CPT violation, the two- 
generation approximation simplifies the expression (39) 
for the first-order oscillation probabilities. With the two 
flavors being and Vr, we have 



Si'J = 5W = 5W = = 0. 



(64) 



This implies no mixing with electron neutrinos occurs in 
the first-order perturbation. Also, inspection of the form 
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of (A^iV)cd(i) given in Eq. (32) reveals that only those 
coefficients for Lorentz violation lying in the two-flavor 
{/ur} subspace can lead to first-order effects. 

Explicitly, we find the first-order oscillation probabili- 
ties are 




= 2LIm((5(°))*wW) 

w Re((5V)L sin(Am32L/2£;). (65) 

Wc assiinie maximal mixing in the last expression, in ac- 
cordance with the parameter values (63). Note that only 
the real part of ((5/i)^x contributes to first-order mix- 
ing in this limit. Also, the corresponding antincutrino 
mixing probabilities are found by the index replacements 
{/xr} — > {/If}, which is equivalent to changing the sign 
of the coefficient (0^)^^^ in {Sh)^r- 

The oscillation probability (65) can be decomposed 
into sidereal amplitudes according to Eq. (48). The am- 
plitudes take the form (49) with the definitions (47). The 
two-flavor approximation makes it straightforward to ex- 
press the latter directly in terms of the real parts of the 
fundamental coefficients (ol)^,- and {cl)'^t rather than 
the intermediate combinations (ol)^!- and (cl)^^. We 
can use these results to estimate experimental sensitivi- 
ties to the real parts of (aL)J^^ and (cz,)^^ for any speci- 
fied experiment. 

Table III presents the results of estimates for the 
eight long-baseline beam experiments considered above. 
In each experiment, numerical values are listed for the 
weighting of the real parts (ol)"^ and (cl)'^^ in the four 

amplitudes {P^^J^.r, (^1'^),.^, {P^]^)^.r, and (P£^)^r. 
The last row provides a rough approximation to the 
attainable sensitivity, based on assuming that a 10% 
sidereal variation in the oscillation probability can be 
detected. The results indicate that all these experi- 
ments can achieve impressive sensitivities to perturbative 
Lorentz and CPT violation. 



C. CP and CPT asymmetries 

In this subsection, we discuss the effects of Lorentz 
violation on tests of the discrete symmetries CP and 
CPT using neutrino oscillations. Experimentally, nature 
is known to break CP invariance in the weak interactions, 
although no CP violation in neutrino oscillations has yet 
been detected. In contrast, compelling evidence for vio- 
lation of CPT symmetry in any system is lacking to date. 
On the theoretical front, CPT invariance has a profound 
connection to Lorentz invariance in quantum field theory, 
where the CPT theorem shows that under mild assump- 
tions CPT violation is accompanied by Lorentz violation 
[8]. No such relationship exists for CP, and indeed CP 
may be violated even when Lorentz invariance holds. 

Consider first CP violation. The CP transformation 



interchanges oscillation probabilities according to 

P,„^,, SPp^^P^, (66) 

so CP violation can be revealed as differences in neu- 
trino and antineutrino probabilities. A generic measure 
of CP violation for mixing involving fiavors {a,b} is the 
asymmetry 

•^ab — p MP--' ^ ' 

which is zero when CP is a symmetry of the oscillations. 

Similar results hold for CPT. Under the CPT trans- 
formation, the probabilities exchange according to 

P.„^.,°^.^Pp,^p„. (68) 

For mixing involving flavors {a, 6}, we can therefore de- 
fine the asymmetry 

= p'^'Tp-^^' - (69) 

Pa— >P!) ' -'P!)— >Po 

This asymmetry vanishes if CPT invariance holds. Note, 
however, that the converse is false: models can be con- 
structed in which CPT is violated even when the asym- 
metry A'^f^ vanishes [13]. In such cases, detailed studies 
of energy and direction dependences may be required to 
reveal CPT violation. 

For the special case of two-flavor models, the probabili- 
ties are blind to possible T violation, so P^,^^^^ = Pi^^^^a 
and Pp^^jy^ = Pp,^^p,^ . Consequently, the CP and CPT 
asymmetries are identical in the two-flavor limit, 

-4Sr=-4Sr^ a,b = ^,,r. (70) 

Note, however, that CPT violation may present itself in 
other ways, such as unconventional energy and direction 
dependence. 

The above asymmetries can be used to test both CP 
and CPT in neutrino-oscillation experiments. In the 
Lorentz-invariant case, violation of CP symmetry occurs 
when both the mixing angle ^13 and the CP phase S are 
nonzero. Measuring ^13 and searching for CP violation 
are major goals of many forthcoming oscillation experi- 
ments. Some experiments can change polarity, choosing 
to focus either positively or negatively charged mesons 
into the decay pipe, and hence can run in both neutrino 
and antineutrino modes. This feature may permit high- 
statistics direct searches for CP violation. The nature of 
the beam or other properties may also lead to accumula- 
tion of neutrino and antineutrino data. In all these cases, 
both Lorentz-conserving and Lorentz- violating situations 
can be accessed, thereby enabling also searches for CPT 
violation. 

The interpretation of asymmetries constructed from 

experimental data requires a theoretical framework. One 
phenomenological approach to CPT violation assumes 
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different masses and mixing angles for neutrinos and an- 
tineutrinos. In the two-flavor case, for example, this 
approach takes a set of parameters {Am'^,6) for neutri- 
nos and a second set (Am^,^) for antineutrinos. It is 
tempting to adopt the resulting explicit expression for 
the asymmetry A^f^ = A'^f for purposes of data anal- 
ysis and interpretation, treating the parameters Am^, 
Am^, 9, and 9 as Lorentz-scalar constants. However, 
according to the CPT theorem this procedure is incon- 
sistent with quantum field theory because under mild as- 
sumptions CPT violation in field theory must come with 
Lorentz violation [8], so the parameters Am^, Am^, 0, 
and 9 cannot be Lorentz scalars. Instead, they must de- 
pend on the 4-momentum of the neutrino, including both 
the energy E and the propagation direction p relative to 
the Sun-centered frame. A typical experiment involves 
neutrinos spanning a spectrum of values for E and p. 
The 4-momentum dependence of the asymmetry there- 
fore entails significant consequences for data analysis and 
its interpretation in searches for CPT violation. 

As an illustration, we derive here the explicit first-order 
form of the two-flavor CPT asymmetries A^^^ = A^^ 
and w4!7,f ^ = A^,f^ in the field-theoretic context. For def- 
initeness we assume maximal mixing, which is consistent 
with the parameter values (63). At first order, calcula- 
tion reveals these asymmetries depend on the coefficients 
(az,)^^ for Lorentz and CPT violation but are indepen- 
dent of (Ci)«^. 

To present the asymmetries, it is convenient to in- 
troduce the CPT-odd part {5h)'^^^ of the perturba- 
tive hamiltonian {Sh)nr with coefficients expressed in the 
Sun-centered frame, 

-{N^'iaL)^^ + {aL)l,) coswe^®. 

(71) 

In terms of this quantity, we find that the CPT asymme- 
try A';^P^ is 

-A^r = - 2i cot (^^)Re {6hf,r- (72) 

This result is valid provided the experiment operates 
away from the region of parameter space leading to small 
oscillations, sin [Aml^L/'iE) ^ 0. For the second CPT 
asymmetry A'^j^'^ , we obtain 

= ^^;«-2itan(^)Re(^/.)^r, (73) 

where now we assume the experiment operates away from 
the region of parameter space leading to large oscilla- 
tions, sin (Am§2i/4-E) 7^ 1. Inspection of these results 
reveals that the two asymmetries A'j^^^ in Eq. (72) and 



A'j^^^ in Eq. (73) contain the same essential information 
about CPT violation but are valid in different regions 
of parameter space. In practice, at least one of the two 
asymmetries can be applied for a given experiment. 

The results (72) and (73) display several interesting 
features. The asymmetries grow with bascilinc L, so ex- 
periments with comparable statistical power but longer 
baselines have improved sensitivity. According to Eq. 
(71), the asymmetries also vary with sidereal time Tg, 
and depend on the direction of the neutrino beam. Both 
these effects are features of CPT violation and its accom- 
panying Lorentz breaking. We remark in passing that the 
structure of the above equations bears a close similarity 
to that of the analogous measures for CPT violation in 
studies of neutral mesons. For example, the dependence 
of {5h)^^'^ on the coefficients (a/,)^^ in Eq. (71) parallels 
that of the measure of CPT violation given in Eq. (14) 
of Ref. [65]. 

In a given experiment, measuring the amplitudes of the 
sidereal variations in the asymmetries (72) and (73) may 
produce interesting sensitivities to the coefficient com- 
binations (iV^(ai);f^ - N^{aL)lr) and {N^{aL)^^ + 
{0'l)^t) ■ These combinations are independent of the 
coeflncients (cl)^^ for CPT -even Lorentz violation. In- 
spection of Eq. (71) reveals that each asymmetry also 
depends on the coefficients (ol)^!- and (a^)^^, which are 
inaccessible via direct sidereal decomposition of the oscil- 
lation probabilities or asymmetries. One way to extract 
sensitivity to these coefficients is to average the data over 
time, in analogy to the extraction of the corresponding 
coefficients for CPT violation in experiments with neu- 
tral mesons [66]. The time-averaged asymmetry A^^'^ 
is 

ACP^ ^2L cot (^^)Re [(a^)^. - iV^(a^)^.] , 

(74) 

while the time-averaged asymmetry A'~l^'^ is 

ACPT „ _2Ltan(^^)Re[(a^)^,-iV^(a^)^,]. 

(75) 

Note that these results remain dependent on the beam 

direction despite the time averaging. Each asymmetry 
therefore typically has distinct physical meanings for dif- 
ferent experiments. For example, the directional fac- 
tor is ~ 0.1 for K2K, ~ 0.6 for MINOS, 
~ -0.4 for OPERA and ICARUS, ~ 0.6 for 
NOz/A, ~ -0.01 for T2K, ~ 0.2 for DUSEL, 
and ~ -0.1 for T2KK. 

While experiments capable of CP tests necessarily test 
for CPT signals in the two-ffavor approximation, the 
CPT signature P^^^^^ ^ Pp^^p^ may be more challeng- 
ing to detect in three- neutrino scenarios. Data from the 
accelerator experiments discussed above or from next- 
generation studies using a beta beam [60] or a dedicated 
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neutrino factory [67] could be well suited for seeking 
three-flavor CPT violation through direct comparisons of 
neutrinos and antineutrinos using the asymmetry A'^^^ 
of Eq. (69). Suitable comparisons of neutrinos and an- 
tineutrinos, perhaps including time averaging as above, 
could also lead to measurements of coefficient combina- 
tions without accompanying sidereal variations. 

IV. COEFFICIENTS g^^ AND 

In this section, we discuss the dominant effects on neu- 
trino oscillations arising from the coefficients '5°'^ and 

H^i^. As shown in Eq. (39), these coefficients leave the 
oscillation probabilities unaffected at first order. The 
dominant effects appear at second order, where the prob- 
abilities are given by Eq. (41). 

The features introduced by q^i and H% include uncon- 

f "ab ab 

ventional energy and directional dependences. However, 
some key differences arise compared to the case of the 
coefficients (0^)"^ and (cl)"^ . For example, the dom- 
inant sidereal variations include higher harmonics with 
frequencies up to 4a;0. Another example is mixing be- 
tween neutrinos and antineutrinos [13], which violates 
lepton-number conservation. This feature arises because 
g^^ and H^j^ lie in the off-diagonal blocks of the pertur- 
bative hamiltonian (9). 

In Sec. IV A, we focus on the second-order contri- 
butions to the neutrino-antineutrino oscillation proba- 
(2) 

bility PpJ-^u^: which involves lepton-number violation. 
The second-order effects on the neutrino-neutrino and 
antineutrino-antineutrino mixing probabilities P^J-^v^ 
and Pp^i^j;^ are considered in Sec. IV B. 

A. Oscillations violating lepton number 

In this subsection, we derive the sidereal behavior 
of the second-order oscillation probability PpJ-^^^ for 
neutrino-antineutrino mixing. Equations (40) and (41) 
specify this probability in terms of the perturbative 
hamiltonian ^/i^s, which itself depends on the coefficients 
'g'^^ and H'^j^ according to Eq. (15). Note that neutrino- 
antineutrino oscillations are independent of the coeffi- 
cients {aL)ab ^^'i i^L)ab ^his order. 

To determine the sidereal decomposition of the off- 
diagonal block (^ft„b of the perturbative hamiltonian, we 
note that H% is an observer vector and hence induces ef- 

ab 

fects at frequency w®, while g^^ is a 2-tensor and hence 
induces effects at and 2^0. The sidereal decomposi- 
tion of therefore takes the form 

= iC)ab + (A)a6 sinwffiTe + (A)ab cosweTe 
+{Bs)ab sin 2uj(sT(s + {Be)ab cos 2u;©re. (76) 



In this expression, the amplitudes (C)„5, {As)abj iAc)ab^ 
{Bs)abJ ctnd {Bc)ab direction-dependent linear combi- 
nations of the coefficients q"i and for Lorentz viola- 

" ab ab 

tion. 

The direction dependence is governed by two vectors, 
the momentum p and the polarization e+ . The momen- 
tum is determined by the beam direction, which varies 
sidereally and is specified at time T© = in the Sun- 
centered frame by the vector (7V^, iV^, N^) given in lo- 
cal spherical coordinates by Eq. (46). Wc denote the 
analogous vector for e+ by ,£^,£^). In the same 
local spherical coordinates, this vector has components 

£^ = -^(cosx(cos0cos^ — «sin^) — sinxsin^), 

£^ = ( cos 6 sin + i cos 4>) , 

v2 

£f = — — ^(sinx(cos0cos^ — i sin0) -|- cosxsin6'), 

(77) 

where we have used the expression (18) for {e+Y ■ 

Some calculation reveals that the sidereal amplitudes 
in Eq. (76) are given as 

{CU = -i.V2[£lH^^-£lErj 

^£lN'E(sT - \gl- - irj) 

{AsU = -iV2[-£lH^^+£^Hl 

^£lEg---£-Eg2 

-£lN^Eg^^+£^N^Eg^/ 
-£lN-Eglf^£lN-Erg-y\^ 

(A)„5 = -ir2{£lBl+£lHl^ 
-£^Eg-^-£lErj 
+£^N'Eg^f+£lN^Eg^,^ 
+£fN''Eg^f + ilN^'Eg^J] , 

{Bs)ab = -^V2[\{£iN- -£lN^)E{sY +gY) 

iBe)ab = -^V2[k{£^N^ + £XN-)EigY +9:f) 
m£-N--£XNy)EgY-gY)].i78) 

This completes the decomposition of Shg^jj in terms of the 
sidereal time T©, the coefficients g"^, H"j^, and the com- 
ponents iV-^, £:l. The analogous decomposition for Sh^^ 
is obtained by taking the hermitian conjugate, following 
Eq. (15). 

At dominant order, the neutrino-antineutrino mixing is 

controlled by the linear combinations H.^'^r and T-l-V given 

ab ab o 

in Eq. (40). The sidereal dependence of (J/i^;, transfers to 
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these combinations, leading to the expansion 

+(4i))„5sina;eTe + {A^^^),icosoj^T^ 

+(B(i))„bSin2a;e7^e + (B^'))„5Cos2a;eT^e, 

(79) 

with a similar expression for H^i^ . The coefficients g"^ 

and H"^ appear in this expansion in linear combinations 
weighted by the complex experiment-dependent factors 
(A^^y)cj and (jM^f,')cd given in Eq. (32). It is convenient 
to introduce the definitions 

cd 
cd 

cd cd 

kb = E(^iV)c-d^"d = E(^iV)c-<i^r- (80) 

cd cd 

In terms of these, the sidereal amphtudes {C^^^)^i^ 

iAi'^)at, {Ai'\-„ {Bi\-„ iBi\i take the same form 
as the corresponding amplitudes in Eq. (78) but with 

q"f and H% replaced with and H^r. For the coefR- 

^ah ab ^ ^ ah ab 

cient combinations 'H'^^ , we can define analogous sidereal 

amplitudes (C(i))5„, {Bi\,, (bP)-,,. 

The forms of these can also be obtained from Eq. (78), by 

first taking the hcrmitian conjugates of the expressions 
on the right-hand side and then replacing g^^ , H^f with 

^ba ' ba 

According to Eq. (41), the combinations Ti^^^ con- 
tribute quadratically to the second-order neutrino- 

(2) 

antineutrino probabilities PpJ_^^^ . This implies that side- 
real variations at frequencies up to 4a;0 are observable. 
Consequently, we expand P^^L^^^ as 

^2 — l"-a6 I 

= iP^'\-b + (Pi'\5sina;eTe + (Pjfj)„5Cosa;eTe 
+ (p£^)„bSin2a;eT'e + (P^^s cos 2a;eTe 
+ (4? )a6 sin BweTe + (P£^ U cos 3u;^T^ 
+ (Pi.^^)„5sin4a;eTe + {P^^^),t,cos4ij^T^. (81) 

Each of the nine amplitudes in this equation is a 
quadratic combination of the coefficients g"^ and H"j^ 
for Lorentz violation. These combinations depend on the 
mass matrix and also vary with the experimental sce- 
nario through the neutrino energy and the direction of 
propagation. 



The explicit forms of the amplitudes in Eq. (81) are 
somewhat lengthy and are omitted here. However, we 
can obtain compact expressions in terms of the ampli- 
tudes for the sidereal decomposition of which are 
defined in Eq. (79). For the harmonics up to 2lu^, some 
calculation yields the results 

{PP)ab = l(C(^))„d^ + i|(4^)).d^ + i|(4^^)a5p 

+ i|(SW),,P + i|(6W)„-,p, 

(p^ju = Re[2(c(i)):,(4')),, + (4i)):,(6(i)),5 

-iA^^^UBi''>)ab], 

[pfXb = Re[2(c(i)):,(4'))„,+(4i)):^(s(i))„-, 

HA'^''):M^)ab\, 
{p^%, = Re[2(c(i)):,(sf))„, + (4')):5(4'))„5], 

i4%b = 2Re[(C«):,(^W),,]-|(4i)),5|2 

+ \W^)ab\^ (82) 

while for the harmonics at Sw® and Auj(^ we obtain 

{p^Xb = Re[(4'^):,(^«),5+(4'^):5(^i'^)a6], 

{P^%, = Re[(4i)):,(^W),,-(4i)):5(B(i))„5], 

(Pi?)),5 = Re[(e«):^(B«)„5], 

(4?)a6 = \{Bi'\if-\{Bi'\i\'. (83) 

The structure of these equations reflects the frequency 
dependence in the sidereal decomposition (79) of H^^jj . 

(2.) (2.) 

For example, the amplitudes (P^ )„(,, {P-f )ab for the 

(2) 

fourth harmonic 4a;0 of the probability Pp^^i^^^ involve 
quadratic products of the amplitudes (M^^)ah7 {Bc^^)^^^, 
for the second harmonic 2lu^ of Ti^^j^ , as expected. 
Comparable expressions for the CP-conjugate transi- 

(2) 

tion probability P^J_^p^ can readily be obtained follow- 
ing the same procedure. The results take the same form 
as Eqs. (81) and (83), but with the index replacement 
{ah} {ab}. 

In the event that neutrino-antineutrino oscillations are 
observed in nature, the sidereal decomposition of the 

(2) 

probability P^J_^^^ and its CP conjugate offers a powerful 
approach to identifying the relevant coefficients and 

H^j^. Each experimental analysis separating the avail- 
able sidereal harmonics would generate eight indepen- 
dent measurements, with multiple experiments able to 
constrain much of the available coefficient space. 

B. Oscillations conserving lepton number 

The analysis in the previous subsection demonstrates 

that the detection oi v v oscillations is a unique sig- 
nal for nonzero coefficients q"^ and H%. However, these 

^ ab ab ' 
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coefficients also contribute at second order to the more 
conventional ^ v and v v mixings. For complete- 
ness, we present the associated equations in this sub- 
section. Effects quadratic in the coefficients (ol)"}, and 
{cL)ab ^^so appear at this order. Inspection of Eq. (41) 
reveals that these contribute independently to the os- 
cillation probabilities, so we set them to zero here for 
simplicity. 

The probabilities for ;/ ^ ;/ and v ^ D mixing are 
affected at second order by Sh^^ through its quadratic 

appearance in the quantities H^i^ and h'^j^ defined in 
Eq. (42). This produces sidereal variations at harmonics 
up to frequency 4cj0. We can therefore decompose 
as the sidereal expansion 

cdef 

+ {Bi^'>)abSm2iVQT(s + {Bf^)abC0s2uj(S)T® 
+ (I?f))absin3we7'e + {'Df^)ab cos3w©Te 
+ (^i2))a6 smAu^T^ + {Tl?^)ab cos4weT^©- (84) 

The nine amplitudes in the above expression can be writ- 
ten as combinations of the experiment-dependent factors 
i^ab )cdef ^.nd thc fivc sidcrcal coefficients for the per- 
turbative hamiltonian Sh^j^ listed in Eq. (78). For the 
amplitudes of the harmonics with frequencies 2wq or less 
in the expansion (84), we find the results 

(C(^))„6 = Y.^M^ab)cd-eA{C)ca{C),s 

cdef 

+ i(6,,).rf(6,,)g/ + i(Se),j(Be)e-/], 

{Af\b - E(-^i6)cJg/[(C)cj(A)e-/ + (A)cd(C)e-/ 

cdef 

-i(A)cj(Sc)e-/ + \{Ao)M-ef 
+ i(B.),j(A)e-/-i(Be),j(A)e-/], 

cdef 

+ \{As),a{Bs)-ef + i(A)cd(^c)e-/ 
+ \{B,),i{As),f + i(Se)ed(^c)g/], 

{Bf\b = E(-^i'b')cde-/[(C)cd(e.)e-/ + (fi.)c<i(C)e-/ 
cdef 

+ |(A)cd-(A)e-/ + \{Ac),i{As)-ef], 

{Bf\b = E(^i6)c&/[(C)cj(Sc)e-/ + (Be)ej(C)e-/ 

cdef 

-i(A)c<i(A)e-/ + ^{AcUiAchf] ■ (85) 

For thc remaining harmonics in the expansion (84) with 
frequencies 3w© and 4a;0, the results for the amplitudes 



are 



cdef 

xi[(A)cd(Sc)e-/ + (^c)cj(S.)e-/ 
+ {Bs)cdiAchf + {Be),Ms)ef], 

cdef 

-{Bs)M)ef + {Bc)cd{Ac)ef], 

cdef 

x^[{BsU{B,hf + {B,U{B,hf], 

cdef 

X i [ - {B.UiBshf + {Bc)M,f] . (86) 
Analogous expressions for the sidereal decomposition of 

(2) 

the quantities n-^^ and thc resulting amplitudes can be 
obtained by substituting barred for unbarred indices and 
vice versa. 

The second-order probability for neutrino-neutrino os- 
cillations inherit the same sidereal-frequency structure. 
Introducing the expansion 



d(2) 
i2 



Re ((5I°Vh1^,0 



(2)^ 



= (Pi'^)ab + (P£^)a6sina;er® + (Pjfj)„6 coswe^e 
+ (P^'^ )ab sin 2ij(sT(B + {P^} )ab cos 2a;eTe 
-t- {P!j^] )ab sin 3a;®Te + (P^^j )ah cos 3a;®Te 
+ (Pg) U sin 4a;eT® + (Pj.^^ )ab cos 4a;eTe , (87) 

we find the nine corresponding amplitudes for the prob- 
ability are given by the equations 



(p(^) 


)ab — 


-Re{{S^:Jnc('\b), 


ip'I] 


Jab — 


-Re((5,SV(4^))„,), 




Jab ~ 


-Re((5(?)*(4^))„,), 


(4? 


)ab ~ 


-Re((5(°))*(B(^))„,), 


(4? 


Jab 


-Re((5(°))*(B(^))„,), 




Jab ~ 


-Re{{s!^^nvi'\b), 


(pS 


Jab ~ 






Jab ~ 


-Re ((5(°))*(^f ))„,), 


(4? 


Jab ~ 





(88) 



The probability for antineutrino-antineutrino oscillations 
can be found from the above equations by replacing all 
indices {ab} with {ab}. 
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The calculations in this subsection demonstrate that 
searches for sidereal variations m v ^ v and D ^ P 
oscillations at the higher frequencies 3u!q and 4a;© can 

offer access to the coefficients q"^ and iJ"- for Lorentz 

ab ab 

violation without the need to study v ^ jy oscillations. 
Moreover, in addition to studios based on the above di- 
rect sidereal decompositions, investigation of the CP and 
CPT asymmetries (67) and (69) introduced in Sec. Ill C 
provides another avenue for data analysis. As before, the 
time-averaged versions of these asymmetries offer sensi- 
tivities to coefficients that are challenging to detect in 
searches for sidereal variations. In all these studies, the 
second-order effects enter in conjunction with a factor of 
L'^, so the large baselines associated with the experiments 
considered in Sec. Ill imply that their intrinsic sensitivi- 
ties to the coefficients g^^ and H^j^ arc only mildly sup- 
pressed relative to the sensitivities to {aL)2b and {cL)ab- 

V. SUMMARY 

In this paper, we study the effects of perturbative 
Lorentz and CPT violation on neutrino oscillations dom- 
inated by mass mixing. The primary focus is on correc- 
tions arising from renormalizable operators for Lorentz 
violation within effective field theory. In the neutrino sec- 
tor, those operators are controlled by SME coefficients for 
Lorentz violation denoted (aL)^^, (cl)^^ , g"^^ , and H^^. 
They can affect conventional oscillations in i/ and 
D ^ D mixing. They also can induce u ■(-^ u mixing, 
which violates lepton number. 

Using time-dependent perturbation theory, a series ex- 
pansion for the oscillation probabilities is derived in Sec. 
II. To second order in coefficients for Lorentz and CPT 
violation, the probabilities for a nondcgenerate mass 
spectrum are presented in Eqs. (35), (36), (39), and (41). 
At first order, only the coefficients (aL)"^ and (cl)^^ 
contribute, and lepton number is preserved. Oscillations 
involving v ^ v mixing appear at second order, governed 
by the coefficients o"^ and H^. 

^ ab ab 

A key feature introduced by Lorentz and CPT vio- 
lation is variations in the oscillation probabilities with 
sidereal time. The sidereal dependence arising from the 
coefficients (a^)"^ and (cl)"^ is discussed in Sec. Ill A. 
It is described by the expansion (48), which involves first 
and second harmonics in the sidereal frequency. At this 
order, the amplitudes for each harmonic are linear com- 
binations of {aL)ab ^'^d i^L)ab ■ Data analyses using bin- 
ning in sidereal time can therefore measure these coeffi- 
cients. 

Section III B addresses the methodology for data anal- 
yses and provides illustrative estimates of numerical 

quantities relevant for sidereal investigations in several 
long-baseline experiments. The results are summarized 
in Tables I through III. For the three-generation case, we 
demonstrate the procedure to identify the relevant lin- 
ear combinations of {aL)2b and {cL)2b ' using the K2K, 



MINOS, OPERA, ICARUS, NOiyA, T2K, DUSEL, and 
T2KK experiments as examples. The two-flavor limit 
is also considered. In this case, the sidereal expansion 
is considerably simplified and the oscillation probability 
takes the comparatively elegant form (65). 

In addition to direct sidereal studies, Lorentz and CPT 
violation can be sought through analysis of CP and CPT 
asymmetries in experimental data. This topic is ad- 
dressed in Sec. IIIC. Suitable CP and CPT asymme- 
tries are defined in Eqs. (67) and (69). In the two-flavor 
limit, these coincide and take the comparatively simple 
form (72) or (73). Experiments running in both neutrino 
and antineutrino modes can probe CP and CPT via this 
route. Analyses along the lines proposed here could pro- 
vide access to different combinations of coefficients for 
Lorentz and CPT violation, including ones that are chal- 
lenging to detect via studies of sidereal variations. 

In Sec. IV, we consider effects arising from nonzero 
coefficients q°'!l' and 11%. Among the features is mixing 

. ab ^ y ^ ^ _ *^ 

between neutrinos and antineutrinos, implying violations 
of lepton number. These coefficients have no first-order 
perturbative effects. Their dominant contributions arise 
at second order, where the probabilities involve quadratic 
combinations oiq"i and H"t. This induces sidereal vari- 

^ab ab 

ations with harmonics up to four times the sidereal fre- 
quency in all three kinds of mixings, u ^ u. ^ v', 
and i/ C". The probabilities for neutrino-antineutrino 
mixing are given in Eq. (81), while those for neutrino- 
neutrino mixing take the similar form (87). 

Overall, we find that the dominant effects of renormal- 
izable operators for Lorentz and CPT violation in the 
neutrino sector generate variations in oscillation prob- 
abilities up to four times the sidereal frequency. Sub- 
dominant perturbative effects may also offer useful infor- 
mation. These higher-order perturbations cause sidereal 
effects at higher harmonics, with signals suppressed com- 
pared to the ones discussed here. We remark that other 
harmonics can also arise from Lorentz- violating operators 
of nonrenormalizable dimensions [13]. A comprehensive 
SME-based study in analogy to that performed for elec- 
trodynamics [68] could establish the corresponding sig- 
nals of Lorentz and CPT violation in neutrinos. 

The results in this work demonstrate that excellent 
sensitivity to Lorentz and CPT violation is attainable 
by studying neutrino oscillations with high energies and 
long baselines. Our primary focus has been beam experi- 
ments, where existing constraints [10-12] span only a few 
percent of the available coefficient space. The procedures 
outlined in this work provide access to essentially all the 
coefficient space, and moreover at sensitivities that can 
exceed the current ones by about two orders of magni- 
tude. 

An interesting direction for further work using a longer 
baseline is a systematic investigation of perturbative ef- 
fects of Lorentz and CPT violation on solar neutrinos, for 
which day-night and annual signals play a role analogous 
to sidereal effects in beam experiments. Future searches 
for Lorentz and CPT violation using extreme baselines 
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could also include studies of oscillations and dispersion 
for supernova neutrinos, for which a sufficient population 
over a substantial solid angle would offer interesting sen- 
sitivity to a significant portion of the coefficient space. 
A more speculative possibility using a cosmological base- 
line would be the search for anisotropies in eventual ob- 
servations of the cosmic neutrino background. The max- 
imal baseline makes this an ideal arena for studying low- 
dimension operators for Lorentz and CPT violation, in 
analogy to the tight limits achieved on low-dimension op- 
erators in the photon sector using observations of the cos- 
mic microwave background [69]. In the meanwhile, the 
long baselines involved in the many current and near- 



future beam experiments on the Earth imply impressive 
potential sensitivities to the effects of Lorentz and CPT 
violation, rivaling the best tests in other sectors of the 
SME. 
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